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ABSTRACT 


In this study, the consolidation of a porous elastic soil stratum in welded contact with elastic medium subjected to 
axisymmetric surface loads has been investigated. The compressible fluid and solid constituents are considered. 
Permeability is different in horizontal as well as vertical directions for porous medium. The Laplace-Hankel transform 
technique is utilized for finding the analytical expressions for the displacements caused by normal circular loading. 
The extended Simpson’s and Schapery’s formula are utilized for the Hankel and Laplace inversion, respectively. For 
numerical computations, the continental, as well as the oceanic layer of the Earth, is considered. The effect of the 
vertical displacement with time for both layers is discussed. It is observed that the compressibility of the solid 
components increases the surface settlement for both layers. It is also observed that for different values of anisotropic 


permeability, the initial and final settlement is the same. 
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INTRODUCTION 


For saturated soil, there is a decrease in value under compressive force, which is known as the compressibility of 
soil. Compression is the process that illustrates the decrease in soil volume under an externally applied load. The 
soil particles squeeze, when stress is applied to a water saturated soil, the water will flow out of the soil. If water 
comes out from pores then air can’t take place. This compression of a saturated soil due to a steady static 
pressure is termed as consolidation. The process of soil consolidation was first discovered by Terzaghi (1). He 
suggested a fundamental approach to the study of a fully saturated soil and expanded the one-dimensional theory 
of soil consolidation. Rendulic (2) generalized this theory to three-dimensional case. Booker and Small (3) 
developed a method for the analysis of the consolidation of a horizontally layered soil using Fourier transform 
method. Booker and Small (4) extended the results of Booker and Small (3) for solving the problem of 


consolidation of horizontally layered soils under both axially symmetric and general surface loading. 


Selvadurai and Yue (5) examined the axisymmetric consolidation response of a poroelastic layer lying 
on a rigid impermeable base caused by circular foundation. Mei et al. (6) considered a cross-anisotropic elastic 
constitutive model for the consolidation of a layered soil using the finite layer procedure. Chen et al. (7) 
presented an exact analytical solution for the axisymmetric consolidation of a transversely isotropic, semi- 
infinite saturated soil. Conte (8) presented an effective technique for the study of coupled consolidation in 
unsaturated soils underneath plane strain with axially symmetric conditions. Singh et al. (9) discussed the quasi- 
static deformation of a half-space subjected to axisymmetric surface loads. Menendez (10) discussed the 


consolidation of an elastic saturated soil with incompressible fluid using the finite element method. Rani et al. 
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(11) studied the influence of negative Poisson’s ratio on the deformation of a half-space for normal strip loading as well as 


axisymmetric normal disc loading. 


Kebli and Merzouk (12) provided an analytical expression of an axisymmetric deformation of an elastic layer 
overlying on a rigid support containing a circular hole due to to uniform loading. Ai and Cheng (13) built up the precise 
integration method for the investigation of consolidation behaviors of multilayered isotropic media. Cheng and Ai (14) 
provided the solution for three-dimensional consolidation of transversely isotropic layered saturated soils in the cylindrical 
coordinate system. Lu et al. (15) developed the reflection-transmission matrix method to study the axisymmetric 
consolidation of a layered transversely isotropic saturated soil. Ai and Wang (16) obtained the solutions for multilayered 


poroelastic soils caused by an axisymmetric effect of Biot's compressibility constraints. 


Ai et al. (17) studied the transient interpretation of a transversely isotropic multilayered half-space subjected to 
vertical loading. Li and Cui (18) presented the solutions for the consolidation problem of saturated soils by pumping. Ai et 
al. (19) investigates the multi-dimensional consolidation of transversely isotropic viscoelastic saturated soils. Closed-form 
solutions for one-dimensional consolidation in saturated soils under external loading has been discussed by Deng et al. 
(20). The consolidation for isotropic viscoelastic media with compressible component because of tangential circular loads 
has been investigated by Ai et al. (21). Wang et al. (22) presented plane strain semi-analytical solutions for consolidation 
of unsaturated soil under time-dependent loading. Lo et al. (23) presented theory for a double-layer system consisting 


upper unsaturated and lower saturated zones caused by surface load. 


The aim of this study is to investigate the consolidation of a homogeneous, isotropic, poroelastic soil stratum in 
welded contact with an elastic medium under axisymmetric surface loading. The fluid and the solid constituents are chosen 
to be compressible. The case of normal circular loading is considered. Laplace-Hankel transform technique is utilized for 
finding the analytical expressions for the displacements. The influence of the vertical displacement with time in the 


poroelastic soil stratum has been discussed. The findings may be employed in Geotechnical and Geophysical Engineering. 
GOVERNING EQUATIONS FOR POROELASTIC SOIL STRATUM 


For axisymmetric consolidation, the governing equations are (24) 


Ou, lou, u, Ou, 1 de ado 
ras =o 2 + = =0 (1) 
or ror fr OZ 1-2v 0r Gor 
OS vem: AO, 1 de ado 
rae ad 7 : ”) 
or r or OZ 1-2v 0z Goz 
where 
Ou, u, OU, Sraete 
e= + + , denotes dilatation. (3) 
Or r Oz 


u,,U, are the displacements in ,Z direction, respectively. denotes the pore-pressure. G denotes the shear 


modulus, @ is the Biot-Willis coefficient, V is the drained Poisson’s ratio. 


The fluid flux g = (4,>o>4,) is given by 
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q, =~%, 00/0r, A = 0, 4, =-H,,00/02 (A) 
where q = (4,>4o>4,) and (%,X:) be the Darcy conductivity in (r,z) direction. 


The constitutive equations are 


V 
On, = 28 | Soo + (—2r) | —ao 
o,, =2Gle e e|-0 


@,, =266., 0, =0;,=0 
¢ =aet aig 
MS) 
where 
2G(v, -v) 


= 2 (i = 2v)\(1 = 2v,) ,is the Biot modulus. (6) 


O,,> O99, O,, are the normal stresses in TI, 0, Z direction, respectively and O,72 90,99 Op, are the 
corresponding shear stresses. c. is the increment of fluid content and V,, is undrained Poisson’s ratio. 
The equation of compatibility is obtained from equations (1)-(3) 
nV’o =GV’e (7) 
where 


(1-2v) 5 _ 8 160 @& 


= V = 
iw 2(1-v)” : ar ror dz" 


is the Laplacian operator. 


The continuity equation is defined as 


0g . 
— =—dw 8 
at q (8) 


The combination of three equations i.e., the continuity equation, the Darcy's law and constitutive equations gives 


the fluid diffusion equation which is written in the form 
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eo lec), oo = Slae+io] 
As Or? ror de oz? ~—s Ot M 


SOLUTION OF THE GOVERNING EQUATIONS 


The n"® order Hankel transform fa (k, Z, s) of f @ Z, s) is defined as 


Tks) ile z,s)J,,(kr)r dr 


ot 8 


so that 


F(r,z,s)= [F(k,z,s)J, (kr)k dk 


pee 


where J, (kr) is the nth order Bessel function. 


The Laplace transform is defined as 


ioe) 


Sy 

— 
= 
N 
N 

>” 
ll 


0 


Application of the Laplace transform to equation (1) gives 


[s+ 10 -7,)a Oa 1 06 ade _ 


2 oe ae, =. 
Or? ror r OZ 1-—2v or Gor — 


-st 
[f(r z,te “dt, s isthe Laplace transform parameter. 
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(9) 


(10) 


(1) 


(12) 


(13) 


where U, » €, © are the variables in Laplace transform domain corresponding to u,, €, O, respectively. 


Taking 1* order Hankel transform of equation (13) which yields 


2- 
a ee aes 


: ; “ka =0 
dz 1-2y G 


(14) 


where U,, €, O are the variables in Hankel transform domain corresponding to U, » &, O respectively. 


Equations (7) and (9) implies 


2 2 
(ast) +e a V’-o=0 
Or ror "Oz Ot 


which gives 


G=Be"+Ce”+Be% +C,e 
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where B, , C, $ B, 5 C, are arbitrary constants. 


Equations (7) and (16) gives 


> 1] —mz mz —kz kz 
eé=—(Be' +Ce™ )+ Be ~+C.e (17) 
a 1 1 ) 3 3 


where B,, C. are arbitrary constants. 


Substituting equations (16) and (17) into equation (14) which gives 
L nk = -k k -k k 
u, = - G(me 2) BP ae +Ce™)+z(B,e % = C.e *)+ Bye 7 +C,e 7 (18) 


Similarly, taking the Laplace transform followed by zeroth order Hankel transform of equation (2) and using 
equations (16) and (17) which implies 


nm 


iu, = Gm? 2) BI = Ce”) + 2(Be” + Ce”) + B,e** —C.e** (19) 


Hence, the complete solution of governing equations for poroelastic soil stratum in Laplace transform domain can 
be written as: 


G(r,zZ,S) = (Be +Ce”+Bie%+ Ce") Jo (kr)k dk (20) 
0) 
é(r,z,S) = |[Z@e~ +Ce™)+ Bye + Ce” |J,(kr)k dk (21) 
} J 
G(nZ6)= x,{(Ber’ +Cje"” + B,e*? + Cees, (kr)k? dk (22) 
0 
G,(r,z,s) = x, |(mB,e"” -mC,e"’ + kB,e” —kC,e) Jy (kr)k dk (23) 
0 
4 = nk B —mz C mz B -kz C kz 
HAT 28)= “Gln 7) ie" +Ce )+2/ oe 7 -Cy,e ) 
0 


(24) 
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Gm? ee - ce™) + Z(B.e™ + Ce) 


(25) 
~ k —MZ MZ es 
salnies)=26| mnt 2 —Ce )—kz(B,e *+C,e™) 

1) 
(3B, -KB, Je [50 -KC, Je" |ai(er dk (26) 


(3B, -KB, Je (56 -KC, Je I (kr)k dk (27) 


where the unknowns B 9 C fi ( J 1D. 4) are calculated by the application of the boundary conditions. 


(e)-3p- eG aE Le) 


B,\_n|,,Q-")s. |(B: 
C,) 2G (v, —v)s |\C, 
Ss, =St(c, —c,)k? =c,(m? -k’) (28) 
SOLUTION OF THE GOVERNING EQUATIONS FOR THE ELASTIC MEDIUM 
The governing equations for the elastic medium are (25) 


2 e 
O'Urn 2 ( Hat) ge 1 Cen _o (29) 
Oz 1—2v, Or 
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2 2 
See ee 1 Cen _g (30) 
or r or OZ 1-2v,, 0Z 
where 
e, =divu = ots aig Mai (31) 
or r OZ 


U,;;,U,,, are displacements in 1, Z direction, respectively. V,, be the drained Poisson’s ratio. 


The constitutive equations are 


OorH = 2Gy Ch rH a 


Over: = 2GyxC rots (32) 


where 


13 


O,,17» Ozzy are the normal stresses in 1, Z direction, respectively and O,,,, is the corresponding shear stress. 


Gy denotes the shear modulus. 


Using equations (29) and (30) and taking zeroth order Hankel transform, a suitable solution is given as 


= -kz 
é,, = De 
H 1 (33) 
where @,, is the variables in Hankel transform domain corresponding to €,, . 
Taking Ist order Hankel transform of equation (29) which yields 
= Zz _ 
tn =|->—>—>— D, + D, |e” (34) 
2(1-2v,,) 
Similarly, zeroth order Hankel transform of equation (30) gives 
= Z - 
i =|-—-__~D, + D, |e" (35) 
2(1- 2v,,) 
and 
1 
{4v,, 3}D, + kD, (36) 


kD = 
*  2(1-2v,,) 


where U,,,, U,,, are the variables in Hankel transform domain corresponding to U,;,, U,,,; respectively. 


D, 5 D, 7 D, are arbitrary constants. 
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Hence, the complete solution of the governing equations for the elastic medium is 


i) Zz 7 
Uy (r,z) =a “20-2,) +D, |e “J, (kr )k dk oD 


\(4v,, -3)+ kz} D, + D, |e*J, (kr)k dk (38) 


Aeon oe e"J,(kr)kdk — (39) 


r 1 
sea ie: e J, (kr)k dk (40) 
0 H 


BOUNDARY CONDITIONS 


Consider a homogeneous, isotropic, porous elastic soil stratum of thickness d lying over a homogeneous, isotropic, elastic 


medium as shown in Figure 1. Assume the upper surface of the porous elastic soil stratum is permeable and the normal 


stress (ox )o is applied on it which implies 
oO =(o De OC, =O, o=O at z=-d (41) 


It is also assumed that along the plane Z = O, the porous elastic soil stratum lying over an elastic medium which 


gives 
C2 FO Oe = Os EAH a Si, at z=0 (42) 
Also, the interface (z = 0) is assumed to be impermeable, we have 
q, =0 at = Z=0 (43) 


Let 


ioe) 


(o,,), = | Noo (kr)k dk. (44) 


0 


Impact Factor (JCC): 8.1695 NAAS Rating: 3.76 


Consolidation of a Porous Elastic Soil Stratum Subjected to Axisymmetric Surface Loading 15 


Porous elastic soll stratum 


Elastic medium 


Figure 1: Geometry of the porous Elastic Soil Stratum Lying Over an Elastic Medium. 


The equations (20), (23)-(27), (37)-(40) together with equations (41)-(43) generate a system of eight equations. 


The constants B i C j ( J] =, 2, 4) ,D ( J =1, 2) can be calculated by solving these equations. After substituting the 
values of B,, C j ( J =|, 2,4) in equations (20), (23)-(27), (37)-(40), the expressions are obtained in the Laplace- 


Hankel transforms domain. The explicit expression for the vertical displacement at (r =Z= 0) is given as 


i, (0,0,s) =- “: |b, | {-2kb, (b, +b, )e* + 2kb, (b, -b,)e™** 
TU 


+2[ b, +(b, —b,.)b, -(b, —B,. + B,, )b, Je" 


J, (Ik) 
12kb.b, |—\—? dk 6 
4 : 2] = (45) 
where 
G'=G,,/G 
D = 8kb;(b, + bg){—b, + (b, — bg )(b, — bs + kb; )}e*" 


+8kb; (b, — bg){—b, +(b, + bs )(b, + b, + kb,)}e™ 
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2 
+|-4 [» — b, (b, + kb;) - “oo | +(b7 — bg) 
m 


2 
|-2%. (b, sf kb, ) ~ 4b, C p.| 7 2(kb; (b, + b,) + 2b,b, e™ 
m 


2 
+|-40 (2 — b, (b, +b) + Eb, |B - b;) 
m 


2 
[22 (b, + kb;) + 4b, C y bs + 2(kb, (b, + b,) + 2b,b, em 
m 


+(b, + b;)(b, + b, —(m-k)b,){2b, -(b, — b,) (by +B,,)iel""" 
+(b, + b,)(b, +b, +(m+k)b, ){2b, -— (b, — by )(b, + b,, jen 
+(b, — b,)(b, +b, +(m+ k)b,){2b, —(b, a b, )(b, 4 b,, tel?" 


+(b, — b,)(b, + b, —(m—k)b;){2b, -(b, + b,)(b, + B,,)}e (46) 


and 
2(1- 1- = 
jo (1-v)(1-v,)}s, . betel! 2v}b, . 
(v,—v)s 2(1-v) 
-4 
i ie NG: Apacer 
2(1-v) 2(1-v) 
2k 2k 
b. = b, =—— 
ome oP ae 
{1 - 2v,,}G’ {2(1-v,)}G' 
b, = 3 — ? 
(3 - 4v,) (3 —4v,,) 


b, = b, +mb, , b,, = b, —mb,, b,, = “(m+k)b,, b, = ~<(m—k)b, 47 
m m 


The integral expression for the displacements is similar to the expression for a poroelastic half-space presented by 


Singh et al. (9) for the case d 00. 
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Numerical Analysis 


In this study, the solution is obtained in the Laplace-Hankel transform domain. The extended Simpson’s formula has been 
utilized for the Hankel inversion. For the evaluation of the Laplace inversion, Schapery’s formula has been used. It is an 
easy and more efficient method for computation proposed by Schapery (26) which has been utilized in the problems related 


to consolidation of soil. It is defined as 
¥(t)=[s¥(s)] ay (48) 


where y (s ) is the Laplace transform of Vit ) and [ | indicates that the value of s is (1/2t). 


s=(1/2t) 


For computation of the vertical displacement at the centre of the disc load, the following dimensionless quantities 


are assumed 


W =(21G/0,)u,(0,0), 


T=(Gz,/P)t, v=0¢,/e,=x,/x, (49) 
Also, the consolidation rate is defined as 
u, (t) 7 u, (0) 
u, (<0) —U, (0) 


For numerical computations, the continental, as well as the oceanic layer, is considered. For the continental layer, 


Consolidation rate = (50) 


the poroelastic soil stratum is chosen to be of Westerly Granite for which v = 0.25, vy,= 0.34, B=0.85 (27) and 
for oceanic layer, the poroelastic soil stratum chosen to be of MHanford Basalt for which 


v =0.30, v= 0.31, B= 0.12 (27). The rigidity ratios for the continental and oceanic layers are G,,/G =2,2 


and G,,/G = 1.76, respectively (28). Assuming the elastic medium is Poissonian i.e. Vy = 0.25. 


Figures 2 and 6 depict the influence of compressibility of solid components on vertical displacement at 
r = Z =O for both continental as well as oceanic layer. @ = 1, denotes the incompressibility of the solid components. It 


is observed that the compressibility of the solid components increases the surface settlement for both layers. 
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log, F 


Figure 2: Influence of the Compressibility of Solid Constituents on the 
Time-Settlement Behavior for Continental Layer. 
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The effect of undrained Poisson’s ratio V,, (-1 S<VSV,S 0.5) on the vertical displacement for the continental 


as well oceanic layer is shown in Figures 3 and 7 for different values of V,. V,, = 0.5, denotes the incompressibility of 


fluid constituent. These figures reveal that for the undrained state (fr > O), incompressibility of fluid constituent 


decreases the vertical displacement but for drained state the vertical displacement has no effect. 


0.36 


0.32 


2 0.28} 


0.24 


0.20 
3 


log, T 


Figure 3: Influence of the Compressibility of Fluid Constituents on Vertical 
Displacement at the Origin for Continental Layer. 
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-4 3 -2 -1 0 1 2 
log, T 


Figure 4: Time-Settlement Behavior for Different Values of Anisotropic 
Permeability Parameter 7 for Continental Layer. 


Time-settlement behavior for different values of anisotropic permeability parameter 7 for the continental as well 


as oceanic layer is represented in Figures 4 and 8. Large value of anisotropic permeability (y ) makes faster the surface 


settlement for both the layers. It is also observed that for different values of anisotropic permeability, the initial and final 


settlement is same. 


0.40 


S 
a) 
n 


Consolidation Rate 
oO 
WwW 
S 


0.25 


0.20 
3 


log, 5T 


Figure 5. Effect of the Compressibility of Fluid Constituents 
on Consolidation Rate for Continental Layer. 
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0.3055 


0.3045 


log, T 


Figure 6. Influence of the Compressibility of Solid Constituents 
on the Time-Settlement Behavior for Oceanic Layer. 
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0.29 
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0.27 
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Figure 7. Influence of the Compressibility of Fluid 


Constituents on Vertical Displacement at the Origin for 
Oceanic Layer. 
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0.3065 


0.3055 


0.3045 


Figure 8.Time-Settlement Behavior for Different Values of 
Anisotropic Permeability Parameter 7 for Oceanic Layer. 


The influence of the compressibility of fluid constituents on the consolidation rate for continental as well as 
oceanic layer is represented in Figures 5 and 9. It shows that the compressibility of fluid constituents have significant effect 
on consolidation rate for both layers. In case of undrained state, the consolidation rate for the incompressible fluid 


constituent is less than the compressible fluid constituent. 


0.34 


0.33 


Consolidation Rate 
So 
Ww 
S 


log,,T 


Figure 9. Effect of the compressibility of fluid constituents 
on consolidation rate for oceanic layer. 


CONCLUSIONS 


¢ The compressibility of the solid components increases the surface settlement. 


¢ Incompressibility of fluid constituent decreases the vertical displacement. 


¢ Large value of anisotropic permeability ( v ) makes faster the surface settlement. 
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For different values of anisotropic permeability, the initial and final settlement is same. 


Compressibility of fluid components has a significant effect on consolidation rate. 
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